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Context. Supercluster A in the extragalactic H ll region NGC 2363 is remarkable for the hypersonic gas seen as faint 
extended broad emission lines with a full- width zero intensity of 7000 km s~^. 

Aims. We explore the possibility that the observed broad profiles are the result of the interaction of a high-velocity 
cluster wind with dense photoionized clumps. 
\ Methods. The geometry considered is that of near static photoionized condensations at the surface of which turbulent 

^ mixing layers arise as a result of the interaction with the hot wind. The approximative treatment of turbulence was 

I ' carried out using the mixing length approach of Cant6& Raga. The code mappings ic was used to derive the mean 

O ' quantities describing the flow and to compute the line emissivities within the turbulent layers. The velocity projection 

in three dimensions of the line sources was carried out analytically. 
C/2 , Results. A fast entraining wind of up to ~ 4300 kms~^ appears to be required to reproduce the faint wings of the broad 

Ha and [O ill] profiles. A slower wind of 3500 km s ^, however, can still reproduce the bulk of the broad component 
and does provide a better fit than an ad hoc Gaussian profile. 
] Conclusions. Radial acceleration in 3D (away from supercluster A) of the emission gas provides a reasonable first-order 

^ , fit to the broad line component. No broad component is predicted for the [Nil] and [Sll] lines, as observed. The wind 

' velocity required is uncomfortably high and alternative processes that would provide comparable constant acceleration 

QQ ' of the emission gas up to 4000 kms~^ might have to be considered. 
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1. Introduction Taylor instabilities, which would be significantly delayed 



in low metallicity gas and in the presence of a very ener- 



O ■ NGC 2363, the largest and most massive H ii region in the g^^j^ Westmoquette et al. (2007a,b,c) observed faint 

^ ; dwarf galaxy NGC 2366 (distance - 3.42 Mpc; Thuan & broad wings in a number of starburst galaxies and suggested 

, Izotov 2005) harbors one of the best documented cases of ^^^^ turbulent mixing layers (hereafter TMLs) on the sur- 

^ . hypersonic gas seen as famt extended broad emission lines f^ce of gas clumps, set up by the impact of the fast-flowing 

H : m an increasing number of giant extragalactic Hii regions ^j^g^g^. ^j^^g (Pittard etal. 2005), might account for this 

>T , (Tenorio-Tagle etal. 1997; Westmoquette etal. 2007a, b,c, phenomenon 

2008). In NGC 2363, this broad component, first reported ^^^-^ ^'^^^^^ quantitatively explore whether sim- 

by Roy etal. (1992) appears as a famt pedestal under the r^^^ ^^^^1^ reproduce the profile shape of the 

narrow Ha H/3 and O ill] lines of the so-called knot A ^^^^ ^. observed in NGC 2363. In order to test our 

(Gonzalez-Delgado et al^l994). This pedestal has an fwhm ^^^^^ ^^^^ ^^^^ ^ ^^^^^ ^^^^^^^ -^^ ^^^^^ 

DUCBR^ 2300 km s (Drissen etal. 2009; hereafter ^^^^^^ of stellar supercluster A. We present succinctly 

''' the observations of the broad faint wings in Ha and [Oiii] 

NGC 2363 is ionized by two massive star clusters la- compare them with simple TML models. A thorough 

beled A (age less than 1 Myr) and B (age 3-4 Myr olcU; analysis of the data set is presented in DUCBR. 
see Drissen etal. 2000). The most intense nebular flux is 
associated with supercluster A, both in terms of "normal" 

emission lines and of the broad emission component. The 2 The observations 

broad component is the focus of this paper. 

The origin of the faint broad emission remains unclear. In order to test our model, we have used spectra of the west- 

Tenorio-Tagle etal. (1997) proposed that it is caused by ern region close to the center of supercluster A described by 

the breakout of fast expanding shells due to Rayleigh- DUCBR, which we briefly describe here. These data were 

obtained with the Gemini Multi-Object Spectrograph's 



^ Knot B, ahhough weaker, contains four Wolf-Rayet stars and Integral Field Unit (GMOS-IFU) attached to the Gemini 
one Luminous Blue Variable (Drissen etal. 2001). North telescope. The two-slit mode was used, covering a 
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Fig. 1. The drawing describing the geometry of the super- 
cluster hot wind within which photoionized condensations 
are distributed. 



5" X 7" field of view, centered on the brightest region of neb- 
ular emission. Two gratings were used: R831, covering the 
6025-6760 A wavelength range with a resolution of 1.4 A, 
and B600, covering the 4090-5400 A wavelength range with 
a resolution of 2.7 A. 

The IFU spectroscopic mode allows mapping of the spa- 
tial extent and geometry of the broad component. Using 
the current data set, DUCBR showed that the intensity of 
the weak broad component in Ha is 2-3% that of the nar- 
row component, in agreement with the results of Gonzalez- 
Delgado etal. (1994). 

3. Modeling turbulent mixing layers 

In the broader context of the interstellar and the intraclus- 
ter medium, Begelman & Fabian (1990) derived a prescrip- 
tion for evaluating the temperature of a turbulent mixing 
layer, which Slavin etal. (1993) used to calculate the emis- 
sion line spectrum of such a layer. Rand (1998) established 
an interesting comparison between their model predictions 
and the spectrum produced by the diffuse gas observed in 
the spiral NGC 891. The models that best reproduce the ob- 
servations required a mixture of TMLs and matter-bounded 
photoionized condensations. More recently, Binette et al. 
(2008; hereafter Paper I) improved earlier models in two 
ways: a) the hot mixing layer and the warm photoionized 



gas are integrated into a single albeit stratified compo- 
nent exposed to an external ionizing source, and b) rather 
than considering a single temperature for the TML given 
by the geometrical mean of the warm and hot phases, a 
temperature structure is derived using the mixing length 
scheme presented in the work of Canto & Raga (1991). The 
possibility of computing not only emission line intensities 
(Binette etal. 1999), but also line profiles, has now been 
implemented to allow a comparison of TMLs with the very 
broad but faint [O ill] and Ha wings observed in NGC 2363. 

3.1. A mixing length approach 

We postulate that both the narrow and the wide emission 
line components of the observed profiles arise from ionized 
gas condensations of low volume filling factor that are mov- 
ing randomly. In a typical extragalactic H ii region with an 
H/3 luminosity of 10*''ergs~"'^, this virial turbulence would 
give rise to a velocity dispersior[3 of ^ 25kms~^ (Melnick 
etal. 2000). In NGC 2363, we postulate the existence of a 
fast hot wind originating from the ionizing cluster, which 
fiows around these condensations, giving rise to a mixing 
layer at their surface, as depicted in Fig.[TJ It is this inter- 
action between the fast wind and the condensations that 
would give rise to the observed broad component. To de- 
scribe the mixing layer, we assume for simplicity an infinite 
plane interface, along which a hot wind of temperature T.^] 
is flowing supersonically with velocity Vw with respect to a 
static warm gas layer of hydrogen number density Uq and 
temperature Tq. The turbulent mixing layer that develops 
between the hot and warm phases has a geometrical thick- 
ness lo, as described in the diagram of Fig. [51 The TML 
structure, the static photoionized gas and the hot wind are 
all isobaric, with a pressure Po (~ 2knoTo). The nebu- 
lar optical emission lines take their origin within the warm 
section of the mixing layer as well as within the static pho- 
toionized layers. 

Within the TML layer, between y = and y — lo, 
the gas is entrained and accelerated. From the theoretical 
point of view, the problem of entrainment in a mixing layer 
involves the description of a turbulent flow. However, fol- 
lowing the work of Hartquist etal. (1986), Canto & Raga 
(1991) developed an approximate treatment of this problem 
based on a mixing length approach of 'turbulent viscosity' 
and on laboratory measurements. Canto & Raga (1991) 
proposed a "single parcel" model, considering mean values 
for the flow variables (density, temperature, velocity) aver- 
aged across the width of the mixing layer. Noriega-Crespo 
etal. (1996) proposed an alternative approach, developing 
a model which resolves the cross section of the mixing layer 
(in this model, the variables correspond to averages along 
the z-axis, aligned with the direction of the shear flow). 
This latter approach was explored further by Binette et al. 
(1999). The adoption of such 'mean' flow characteristics al- 
lows us to get around having to deal with the details of the 
turbulence cascade. 

For the case of a thin, steady state, high Mach number 
radiative mixing layer, the advective terms along the direc- 
tion of the mean flow can be neglected with respect to the 
corresponding terms across the thickness of the mixing layer 

^ It is a well documented fact that the line profiles of giant H ll 
regions greatly exceed the width given by thermal broadening 
alone (c.f. Melnick etal. 1988, 2000 and references therein). 
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Schematic diagram showing the cross-section along the y-axis of a plane-parallel mixing layer of thickness lo- The 
layer is formed by the interaction of a supersonic wind (velocity and temperature T^) with a static layer of 
as of H number density Uo and temperature To (at y = 0). 



(Noriega-Crespo et al. 1996, Binette et al. 1999). Under this 
approximation, the momentum and energy equations can 
be written as: 



= 



" dy^ 



dv 



= L-G 



(1) 



(2) 



where ?/ is a coordinate measured from the onset of the 
mixing layer (see Fig. [5]), v the bulk flow velocity (the ve- 
locity component projected along the x axis) as a function 
oi y , L and G are the radiative energy loss and gain per 
unit volume, respectively, and /i and k are the turbulent 
viscosity and conductivity, respectively, which are assumed 
to be constant throughout the cross-section of the mixing 
layer. Eqn.[T]can be integrated to obtain the linear Couette 
flow solution: 



(3) 



where lo is the thickness of the TML and is the velocity 
of the hot wind, parallel to the static layer. As the gas enters 



the TML and starts flowing, its 'mean' velocity in the x 
direction increases linearly until it reaches the wind velocity 
at y — lo- This solution can be substituted in Eqn.[21 which 
can then be integrated to derive the temperature behavior 
across the layer. 

Within the TML, we consider the equations governing 
the fractional abundance /,; of each species i. This abun- 
dance must satisfy the equation: 



D 



dy^ 



(4) 



where 5^- is the net sink term (including coUisional ion- 
ization, radiative and dielectronic recombination, charge 
transfer, etc., and the reverse processes which populate 
the current species) of the species i. The turbulent diffu- 
sivity, D, is of order unity and assumed to be position- 
independent. At the inner and outer boundaries of the mix- 
ing layer, the ionization fractions are set by the equilibrium 
values. 

To complete the description of the mixing layer, we re- 
quire lateral pressure equilibrium (which determines the 
density of the flow along y), and calculate the turbulent 
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viscosity with a simple, mixing length parametrization of 
the form: 



The impact of photoionization can be inferred from the 
behavior of F across the layer, where F is defined as follows: 



H = a pslo , 



(5) 



where p and s are the mass density (gcm"'^) and sound 
speed (respectively) averaged over the cross-section lo of 
the mixing layer (Fig. [2]). From the work of Canto & 
Raga (1991), the value of the proportionality constant a 
is 0.00247, as described in Appendix A of Paper L It is the 
required value for a supersonic mixing layer model to match 
the opening angle of 11° of subsonic, high Reynolds num- 
ber laboratory mixing layers in the limit in which the jet 
Mach number tends to one. 

Considering that the turbulent conduction and diffu- 
sion Prandtl numbers are of order one, we can compute 
the conduction coefficient (where Cp is the heat 

capacity per unit mass averaged across the mixing layer 
cross-section) and the diffusion coefficient as 13 ~ p/'p. In 
this way, we obtain a closed set of second order differential 
equations ([5] and 3]), which can be integrated with a simple, 
successive overrelaxation numerical scheme. 



3.2. TML calculations with the multipurpose code 

MAPPINGS IC 

We use the code mappings ic (Ferruit et al. 1997) to com- 
pute the radiative energy loss term L and the photoheating 
term G (Eqn.[2]) at each position across the TML. At both 
the inner and outer boundaries of the mixing layer, we as- 
sume equilibrium ionization of the different species, while 
across the layer, our simple overrelaxation scheme is used 
to determine the ionization fractions (Eqn.d]). For the ion 
diffusion of each species fi, the spatial differential equa- 
tions are converted to temporal equations, with the use of 
pseudo-time steps At — Ay'^ lo/ (a n s), where n is the aver- 
age H number density. This allows us to use the temporal al- 
gorithm previously described in Binette & Robinson (1987) 
for determining the spatial diffusion of the ionic species. 

The radiative transfer is determined by integrating 
(from the hot layer y — lo down to y = 0) the intensity 
of the UV diffuse field produced by the layer, assuming the 
outward only approximation. Any UV radiation impinging 
the layer from the outer boundary is simply added to the 
diffuse field at the onset of the integration. The intensity 
of the external ionizing field is defined by the ionization 
parameter as follows: 



Uo = ifw/cUo , 



(6) 



where c is the speed of light, ip^ the flux of ionizing photons 
impinging on the TML, and rio, the H density at y = 0. The 
flux ipu] is the sum of two components: resulting from 
external UV sources such as hot stars and ipx due to the X- 
ray diffuse field generated by the hot wind (of column Nx 
and average temperature Tx). In circumstances where the 
hot wind occupies a large volume, the soft X-ray luminosity 
of the hot wind can become significant despite its density, 
which is much lower than the optical line emitting gas. In 
this paper, we neglected the term ipx because the hard 
radiation appears to play a negligible role in comparison 
with the UV from the hot stars. In fact, using the diag- 
nostic line ratio diagrams of Veilleux & Osterbrock (1987), 
we find that NGC 2363 occupies the same locus as that of 
H ii-region like galaxies (DUCBR, Gonzalez-Delgado et al. 
1994). 



L-G 
L + G ■ 



(7) 



The quantity F is zero when the temperature corresponds 
to the equilibrium value and it is unity when cooling dom- 
inates. In our calculations, we iteratively determine To at 
y = until the condition F = is satisfied. 

For the calculations presented in this work, the mix- 
ing layers occur at the outer surface of gas condensations 
immersed in the radiation field of hot stars permeating su- 
percluster A of NGC 2363. The spectral energy distribution 
(sed) that we adopted for ip^ was calculated using the code 
LavalSB (Dionne & Robert 2006). It corresponds to a newly 
formed star cluster 1 Myr old, in agreement with the upper 
age limit derived by Drissen et al. (2000). The stellar masses 
are represented by a Salpeter distribution with an upper 
mass cut-off of 100 Mq. The abundances of the atomic el- 
ements of the gas are set at 20% solar, in line with the 
conclusions reached by Luridiang0 etal. (1999). 

After computing the emission line spectrum of a given 
TML, MAPPINGS IC offers the option of computing sepa- 
rately the emission lines generated by the inner photoion- 
ized layer (i.e. y < in Fig. [5]) where equilibrium ionization 
prevails. A simple isobaric photoionization model is cal- 
culated in this case using the radiation field that has not 
been absorbed by the mixing layer. The total line spectrum 
is then given by taking the sum of the line intensities from 
the TML model (the broad profile component) and from 
the static photoionization model (the narrow profile com- 
ponent). 

Since the mixing layer is isobaric, its density profile as 
a function of thickness is uniquely determined by the pres- 
sure Pq. Given the value of Uo at the boundary y = 0, 
the pressure is derived from To, which is the equilib- 
rium temperature in the photoionized case. Within our se- 
lected 1.2" X 1.2" aperture (Sect.glT]), the observed [Sill 
(6716/6731) ratio is 1.36, which translates into a densitjO 
of 89 cm"'^, assuming a temperature of 12 000°K. Within 
a factor of two, this density is also representative of the 
layers that emit [Nil] or [On]. Thus all models presented 
in this work were calculated assuming no=100cm~'^. It is 
much lower than the critical density of most atomic tran- 
sitions typical of TMLs, which means that the gas density 
perse is not a significant parameter in these calculations. 
In this case, equal external parameters can be considered 
equivalent whenever the product of the H density Ug and 
thickness lo is the same. Therefore, it is sufficient to specify 
the quantity — nolo to uniquely define a model, when 
Uo is kept constant. 

To summarize, in order to compute solutions to the mix- 
ing layer, we must specify the values of the following param- 
eters: the ionization parameter, the mixing layer's nominal 



^ The abundances inferred by these authors are 25% solar for 
the stellar atmospheres and 20% for the nebular gas. 

* This density only applies to the photoionized layer of the 
condensations. The core of the condensations are expected to 
be cold and therefore much denser. 
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columr0 No and finally the temperature and velocity of the 
hot wind: and Vw 

3.3. A high-velocity hot stellar wind 

The full-width at zero intensity (fwzi) of the Ha and 
[Oiii] profiles reaches the remarkable value of 7000 km s~^ 
in NGC 2363. This means that for a geometry in which the 
emitting gas moves radially in 3D, the velocity must ex- 
tend at least up to 3500 km s~^. In practice, faster winds 
are required in our models, because the layers that have 
velocities approaching that of the wind do not produce 
any optical lines. The main reason is that these layers are 
very hot and their densities so low that their line emissiv- 
ities become negligible. Another reason is that these gas 
layers are overionized. Towards the static layers, the tem- 
perature is close to being isothermal because photoheat- 
ing equals radiative cooling. The transition between the 
isothermal layers and the wind dominated layers is quite 
abrupt. Therefore the thermal structure of TMLs consists 
of two zones: the isothermal warm zone and the hot thermal 
bump. Examples of such a thermal structure is illustrated 
in Fig. [313 as a function of the normalized thickness y/lo- 
The three models shown are described in detail in Sect. 13. 41 
below and differ only by their value of No (hereafter N2Q, 
in units of lO'^'^cm^^). The normalized thickness maps di- 
rectly into the velocity domain, since the Couette flow so- 
lution implies a linear velocity increase between the static 
layer and the wind flow (Eqn.[3]). If we define e as the frac- 
tion of the layers' thickness where T^, is isothermal, which 
turns out to be where the optical lines are produced, it fol- 
lows that a wind velocity of 3500 e^^ kms~^ is required to 
cover the velocity span observed in NGC 2363. We could 
not get credible calculations that had e > 0.9. The model 
with No = 33 in Fig.[2l3 has e = 0.81. Hence, in order to re- 
produce the observed broad wings, it is necessary to adopt 
a wind velocity as high a,s Vw — 4300 km s~^. This value is 
somewhat higher than the 3500 km s~^ value calculated by 
Sternberg etal. (2003) for 03 stars. 

Rather than a single stellar wind, we propose that we 
have a wind from a dense cluster of massive stars. As shown 
by Canto et al. (2000), the resulting cluster wind has a ter- 
minal velocity equal to the velocity of the winds from the 
cluster stars. Also, the interaction of the different winds re- 
sults in a very high initial temperature for the cluster wind 
(~ 10^ °K for the wind velocity we are proposing), but this 
temperature rapidly drops beyond the outer radius of the 
cluster, as the cluster wind approaches its terminal velocity. 

Another alternative to a stellar origin for the wind is 
that it has been generated by supernovae. There is, how- 
ever, no evidence of past supernovae in knot A although we 
cannot entirely rule it out. In knotB, there is evidence of a 
cavity surrounding the central stellar cluster, which might 
have been caused by supernovae explosions. The explana- 
tion provided by Drissen etal. (2000, 2001) is that super- 
cluster A is extremely young (^ 1 Myr) while supercluster B 
is older (~ 3-4 Myr). 

We found that the temperature T^, of the hot wind 
is not a critical parameter of the TML models. Its value 



^ The quantity A'o {~ riolo) is a convenient model descriptor. 
However, it is a bad estimator of the true integrated H column, 
which is a lot smaller since the density is not constant but de- 
creases as the temperature rises with thickness y (Fig.|3}. 
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Fig. 3. Behavior in panel a of the temperature normalized 
to the value in the wind and in panel 6 of F (Eqn.[7l), both 
as a function of the normalized thickness yjlo- The three 
models shown vary only by their mixing layer nominal col- 
umn: 7V2°o = 4, 10 and 33. 



can be varied and still result in a sequence of equivalent 
models (i.e. with similar line ratios), provided remains 
the same. However, to ensure that the X-ray flux from the 
hot wind would not be excessive, we have set its tempera- 
ture to Tu, = 5 X 10^ °K in all our calculations. The non- 
detection of NGC 2363 by the satellite XMM-Newton in 
the 0.1-lOkeV X-ray band has allowed the determination 
of an upper limit of 4 x 10~^^ ergcm"^ s~^ (Y. Krongold & 
E. Jimenez-Bailon, private communication). Using the code 
MEKAL (Mewe etal. 1986; Liedahl etal. 1995) and assum- 
ing a cylinder of hot gas of radius 10 pc, length 100 pc and 
density n^, = 100 cm-^ x 12 000°K/T„ = 2.4 cm'^, one of 
us (Y. K.) recently determined that the wind would emit a 
flux more than an order of magnitude lower than the above 
upper limit. Thus, even much higher wind temperatures 
(^ 10* °K) would neither produce a flux that exceeded this 
upper limit. 

3.4. Internal structure of three test models 

After deflning the stellar SED, density and wind character- 
istics, we can now proceed with the calculations. In order 
that the line ratios of the photoionized layers could match 
the nebular values, we selected a high ionization param- 
eter Uo = 0.02. The total [Oiii]/H/3 ratio, including the 
static gas emission, is « 8.1, which is comparable to the 
observed ratio of 8.9. Our main aim is to reproduce the 
broad lines profiles and therefore no further analysis of the 
line ratios of the cores of the profiles has been attempted. 
In Fig. [31 apart from we also show the behavior of the 
measured imbalance between cooling and heating, F (see 
Eqn.[7]), as a function of normalized thickness for each of 
the three models whose mixing layer's column differ and 
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take on values of A^2o ~ 4' ^^"^ ™ units of 10^° cm~^. 
For each model, the ionization parameter is Uo = 0.02 and 
the wind temperature and velocity are T^, = 5 x 10^ °K and 
— 4300 km s^^, respectively. The equilibrium tempera- 
tures at the onset of the layer are To — 11 460, 11 330 and 
12 370 °K for the models with N^q = 4, 10 and 33, respec- 
tively. Because e increases with A^joi these TML models re- 
sult in a broad profile whose width increases monotonically 
with N2Q. Alternatively, a sequence of increasing profile's 
width can be obtained by increasing Vw while keeping 
constant. 

It is well known that the cooling rate L of a hot plasma 
decreases with temperature when its value lies above « 
2 X 10^ °K (e.g. Fig. 8 in Ferruit etal. 1997). This behav- 
ior of the cooling is responsible for the temperature bump 
towards the interface with the wind, which is apparent in 
Fig.[3l In effect, near the wind interface, the heating from 
turbulent dissipation overwhelms radiative cooling and re- 
sults in r near unity. It also results in an overionized plasma 
of very low density (hence low emissivity), which does not 
generate significant line emission, at least in the optical 
domain. The impact of heating by turbulent dissipation in- 
creases considerably when we consider thinner mixing lay- 
ers, as can be seen in Fig. [3] Towards the colder gas lay- 
ers, radiative cooling becomes sufficiently strong to enable 
a balance between radiative cooling and heating by pho- 
toionization and turbulent dissipation. Although heating 
by photoionization dominates in the near-isothermal region, 
turbulent dissipation is nevertheless present as indicated by 
r, which reaches 0.3 in the case of the = 4 model. 
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[Oiii]A5007 (dashed line) as a function of radial velocity 
is plotted for each of the three models of Fig. [31 Note that 
the [Oiii] fluxes have been scaled down by a factor 5. The 
He II A4686 emission flux is shown in the case of the A''2o=33 
model only. 



3.5. Projection in 3D and line profile calculations 

In order to derive the integrated line profiles from the 3D 
distributions of gas condensations, we first compute the 
profile from a single condensation at the surface of which 
the fast wind generates a TML. We neglect the details of 
the contours of such condensation and consider that the 
TML takes place on the sides of the condensation, as if it 
was a cylinder whose axis is oriented radially with respect 
to the wind source. In Fig.[3]we show the resulting emission 
flux from the same three models of Sect. 13. 41 as a function of 
radial velocity. The vertical scale is arbitrary between mod- 
els, but remains the same for lines of the same model. These 
ID profiles correspond to a single condensation that lies ex- 
actly along our line-of-sight to the wind source. The profiles 
have been initially smoothed in mappings ic by a narrow 
gaussian of 25kms^^ dispersion (see Sect. 13. ip . which ac- 
counts for the virial broadening due to the supersonic mo- 
tions of the condensations. To a first order and for most 
optical lines, the profiles resemble a top-hat function that 
extends up to a limit, Vum. The velocity limits up to which 
the Ha or [O ill] flux extends are Vum — —1900, -2800 and 
-3440kms-i for the models with N^q = 4, 10 and 33, re- 
spectively. This top-hat shape is the result of the Couette 
flow solution, which corresponds to a velocity that increases 
linearly with y (Eqn.[3]). The Hen line is an exception be- 
cause it is mainly produced at the onset of the thermal 
bump. Viim is lower than since the optical line emissivi- 
ties become negligible inside the overionized thermal bump 
region, as discussed above. Essentially, all of the flux origi- 
nates from the near- isothermal region (depicted in Fig. [3]). 
The thinner the TMLs gets, the smaller becomes the region 
where emission takes place, which results in proportionally 



narrower line profiles. Furthermore, the model with A'o=33 
is ionization-bounded, while the other two with A^o=10 and 
4 are matter-bounded, which result in significantly higher 
[Oiii]/H/3 line ratios, of 10.8 and 15.7, respectively. 

The emission profiles of Fig. [5] apply to the simple one- 
dimension case. We will now consider the more realistic case 
of a radial wind in 3D and compute the profile that would 
be seen by a distant observer. To achieve this, we assume 
some simplifications about the problem, namely, that all 
the TMLs are immersed in a radial constant velocity wind, 
the wind is isotropic with respect to the wind source and 
the covering factor C/ of the mixing layers is also isotropic. 
To simplify things further, we will assume that the TMLs 
are located at the same radius R from the central source 
of the wind. For any emission line i at rest-wavelength A^, 
these considerations lead to the following integral for the 
fiux: 

Fi=2T:R^Cf< F\v')smede 

+ \ F\-~v') exp(-TA. cos 6*) ^x^ddQ L 

° J (8) 

F\v') = sec 6* FV cos 6*) , 
v' =v cos d , 

z; = c(A - K)/Xl , 

where c is the speed of light, the angle with respect to the 
line-of-sight to the observer, R the radius where the TMLs 
are located with respect to the wind source, and F\v) in 
the integrand is the emission flux for line i as computed 
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by MAPPINGS ic (i.e. for 9 = 0). To include the possibil- 
ity that dust might fill the volume of the sphere of radius 
R and thereby absorb line emission from the (red-shifted) 
TMLs on the far-side, we introduce in Eqn.[5] a dust ex- 
tinction cross-section ct^ evaluated at AJ,. The parameter rid 
is the dust grain density filling the sphere and T\i is the 
resulting line extinction due to dust. The effect of internal 
dust would be to skew the profile towards the blue. Since 
we found no evidence of profile asymmetry, we have set 
TTrf = 0. There could be additional 'intervening' dust that 
covers the whole nebula, however. This possibility is not rel- 
evant to our profile study since the differential reddening 
across the profile width turns out to be negligible, assuming 
a C(H/3) of 0.2 for knot A (Gonzalez-Delgado etal. 1994). 
As for the two limits of integration, whenever the TML sys- 
tem (of projected area ttR^ on the sky) is fully contained 
within the aperture of the spectrograph, these reduce to 
9i = and 6f = 7r/2, which is what is assumed hereafter. 
We have explored relaxing this assumption. With 9f ^ tt/2, 
the profiles become flat-topped while for 6i ^ 0, the profiles 
does not change much in shape but becomes progressively 
narrower. In a more realistic description, the TMLs would 
cover a range in radii and C/ would be a function of ra- 
dius. The current description, however, suffices to capture 
the basic implication of a 3D geometry on the line profiles. 

4. Comparison of TML profiles with the 
observations 

4.1. Characteristics of the broad profile wings 

We analyzed a subset of the available data as follows. To 
derive the highest S/N possible, we considered a square re- 
gion of 1.2" X 1.2" (20 X 20 pc) centered on knot A and ex- 
tracted a red and a blue spatially-collapsed spectra within 
this area. We modeled the underlying continuum and sub- 
tracted it from both the Ha and [O ill] A5007 lines. In order 
to facilitate the comparison of the models with the faint 
broad wings observed in [O ill] and Ha, we converted the 
profiles into velocity space relative to the centroid of each 
line. Direct profile comparison is achieved by simply super- 
posing different lines. 

At the time of writing, various issues concerning the 
absolute calibration of the blue spectrum could not be sat- 
isfactorily resolved. In what follows, we will not rely on 
the absolute flux scale of the blue spectrum, but focus in- 
stead on profile shape comparisons or on relative ratios of 
the broad component with respect to the central core com- 
ponent. After rescaling the blue spectrum until the broad 
[O III] superimposes the broad Ha profile, one finds that the 
'shape' of the broad component is the same in both lines. 
This is apparent in Fig.O where the Ha profile as a func- 
tion of Doppler velocity is shown in red and [O iii]A5007 in 
blue. To express fiuxes, the quantity F^^ = 10^^F„ is used 
hereafter. The Ha line peaks at a value of Fi^=853 (but not 
[O III] which has been rescaled in this figure) . No reddening 
correction has been applied. We now discuss plausible in- 
terpretations of the scaling factors that are obtained from 
superimposing line profiles. 

Detailed analysis showed that if we scaled the [O ill] line 
profile so that the peak of its narrow core equalled that of 
Ha, the broad wings of Ha were brighter than [Oiii] by 
about 50%. If instead of [Oiii], we overlaid the continuum 
subtracted H/? line and rescaled it so that the peak of its 




-4000 -3000 -2000 -1000 1000 2000 3000 4000 

V (km s"') 

Fig. 5. Superposition of the emission line profiles of Ha 
(red fine) and of [Oiii]A5007 (blue line) as a function of 
Doppler velocity with respect to the centroid of both lines. 
The scale applies to Ha and the units are ergcm^'^ km^^. 
The [O III] profile has been scaled up so that its broad wings 
coincided as much as possible with those of Ha. Emission 
lines involving other energy levels have been partially blan- 
keted with dotted line contours. Using a thick gray line, we 
overlaid the 3D profile from a TML model with iV2Q=10 
and T4j=4300kms~^. The yellow line shows how a broad 
Gaussian of 2500 km s""'^ fwhm would compare. The green 
colored dashed-line is a FWHM=110kms~^ Gaussian fit of 
the narrow core component. 



narrow core equalled that of Ha, we similarly found that 
the broad Ha lies w 50% higher than the broad H/3. We 
verified that saturation is not taking place. Flux spilling 
over nearby pixels appears also to be ruled out. A possible 
interpretation is that the gas responsible for the faint broad 
wings is further absorbed by dust than the gas responsible 
for the bright narrow cor^B- Apart from our basic conclusion 
that the broad [O III] and Ha profiles are of similar shapes, 
which suggests that the broadening mechanism acts uni- 
formly on both lines, we find that the data is consistent 
with the broad [O ill] /H/3 fiux ratio being the same as that 
of the narrow core [O ill] /H/3. This follows from the two sep- 
arate comparisons above that indicated that the broad H/3 
and [O III] are depressed by the same factor with respect to 
the broad Ha. It is noteworthy that these two conclusions 
are not affected by our calibration uncertainties. 

Hereafter we will concentrate on the faint broad profiles 
shown in Fig.[5l As discussed above, the [Oiii], Ha and H/3 



® The amount of dust encountered by Gonzalez-Delgado et al. 
(1994) for the narrow line emitting gas in knot A is as little as 
C(H/3)=0.2 (or tv = 0.22). Uncertainties in our absolute cal- 
ibration of the blue spectrum does not allow us to verify this 
value. To have more dust covering the broad line emitting gas, 
however, is somewhat counterintuitive and this particular inter- 
pretation should be considered tentative. 
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Fig. 6. Same notation as in Fig. [5] except for the TML 
model, which corresponds to N2q=3S. For comparison, a 
Gaussian of 3200 km s^^ fwhm is overlayed. 



broad profiles are characterized by a similar functional de- 
pendence on velocity. This similarity in shape and the fact 
that the integrated [O ill] /H/3 flux ratio for the wings is the 
same as for the narrow core imply stringent constraint for 
any proposed model. It suggests that the physical condi- 
tions pertaining to the gas that produces the faint broad 
wings are very likely similar to those of the more quiescent 
nebular gas. 

The FWHM of the core of the profiles was derived by 
fitting a single Gaussian. We obtained 78 and 115kms~^ 
for Ha and [O ill] , respectively. The [O ill] core profile is 
unresolved and appears broader due to the lower resolution 
of the blue grating. In order to distinguish at which velocity 
the broad line flux exceeds the intense flux from the narrow 
core, we overlay in Fig. [5] a Gaussian fit (green dashed line) 
to the narrow core (with FWHM=:110kms~^). Wc infer that 
the observed profiles becomes noticeably wider than the 
narrow Gaussian below « 2, which lies at only 0.23% 
of the peak flux value. Note that the width of the narrow 
Gaussian is 325kms~^ at Fi^=2. 

4.2. Comparison with 3D TML line profiles 

Since our main aim is to reproduce the weak broad wings 
underneath each line, we now focus on the profiles gener- 
ated by the TML, leaving out the narrow line component 
produced by the static layer at the surface of the conden- 
sations. Since the emission flux from the broad component 
is proportional to the surface area of the TML gas exposed 
to the ionizing radiation, this area has to be quite small 
(« 3%) relative to that of the nebular gas (the 'static' or 
narrow line component) . 

We applied the profile integral described in Eqn.[51 
which considers the effects of velocity projection in 3D, to 
each of the three TML models shown in Fig. 01 Because the 
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Fig. 7. Same notation as in Fig.O except for the TML 
model, which corresponds to 7V2o=4. For comparison, a 
Gaussian of 1500 km s^^ fwhm is overlayed. 

narrow line cores in NGC 2363 have a width that exceed the 
value assumed by mappings ic (25kms^^, see Sec. l3.ip . we 
further convolvecQ the calculated profiles using a narrow 
Gaussian of 107 km s"-^ fwhm. The fwhm describing the 
velocity dispersion of the condensations (Fig.[T]) also applies 
to the TMLs. 

The TML model of the [Oiii] line assuming N^q = 10 
is superimposed to the data in Fig. [5] (gray line). It is un- 
necessary to superimpose the Ha model separately, since 
it is visually undistinguishable from the [O ill] model, even 
if the difference in spectral resolution were considered in 
the modeling. The profile fit that is achieved of the broad 
wings, although clearly imperfect, is certainly encouraging, 
given the approximations made to the transfer and to the 
geometry of the TMLs. Because Gaussians are a natural 
profile for describing emission lines over a wide variety of 
physical situations, we also compared our data with sim- 
ple Gaussian profiles. In Fig. [51 the yellow line represents 
a Gaussian with an fwhm of 2500 km s~^. We find that 
our TML model provides a better description of the data. 
(Gaussians with different fwhm are drawn in the subse- 
quent Figures [6] and [7]). There is definitely more flux present 
at both low and high velocities in the observed profiles than 
provided by any single broad 'Gaussian'. The same appears 
to be the case for the TML model with A^o=10, but to a 
much lesser extent. For instance, the TML profile fills the 
widening of the profile beyond the narrow Gaussian core 
(dashed green line) much better compared to the round 
top of a wide Gaussian. A Lorentzian profile would pro- 
vide a better description of the data than a Gaussian, as 
shown in Fig.[8l However, given the physical characteristics 
of the gas, there is no physical justification for a Lorentzian 
profile. 

^ The target width is 115k ms~^, whic h implies a convolution 
by a Gaussian of fwhm of Vll5^ — 25^ = 107kms"\ 



Luc Binette et al.: Broad line wings in nebula NGC2363 



9 



In Figs.[S]and[71 we show the other two TML models of 
equal Vw, but with different of 33 and 4, respectively. 
The model that fits the wings better (Fig. [6]) does worse 
with the core region, and conversely with the other model. 
For comparison, Gaussians that approximately match the 
same region as the corresponding TML model are overlaid. 

To obtain a much improved fit to the profile wings, 
the juxtaposition of two or more TML profiles of different 
widths would therefore be necessary. Because the [O ill] /H/3 
ratio varies markedly along the No sequence (Sect. 13. 5|) and 
since the observed broad [O ill] /H/3 line ratio is the same as 
the one found for the narrow core, we would argue against 
combining models of different thickness No- An alternative 
way to combine TML profiles would be to consider a radial 
gradient in wind velocity Vw Such a sequence in which 
varied and No remained at the value set by the ionization- 
bounded case would present a clear advantage, since the 
resulting [O ill] /H/3 ratio would not change from the value 
characterizing the static nebula case. 

To illustrate this, we looked for an ionization-bounded 
model that would be equivalent to the model with A^jo—l^ 
(Fig. [51 which is the model that fits relatively well the bulk 
of the broad profile). In this new model, the increase in 
thickness is compensated by a reduction in wind velocity. 
The resulting profile is shown in Fig. [5] and corresponds 
to a slower wind of T4,=3500kms~^ and a larger thick- 
ness iV2o=31. The comparison of Fig. [5] with Fig. [S] shows 
that the two profiles are comparable. In the case of the 
slower wind model, the broad profile is slightly wider be- 
cause e is larger (0.88) and it fares better in the wings than 
the 4300kms"^ model with N^q=10. This new model is 
preferable, since it results in the same [Oiii]/H/3 ratio as 
that of the narrow core. Furthermore, the wind velocity is 
more conservative and corresponds to the value calculated 
by Sternberg et al. (2003) for 03 stars. 



4.3. Other emission lines 

The TML models predict that there should not be any 
broad wings underneath the lower excitation lines such as 
[N 11] or [S 11] . The reason is that in models with a high ion- 
ization parameter (e.g. Uo ^ 0.01) the degree of ionization 
throughout the turbulent layer is high and the ionized gas 
is dominated by high excitation species, such as 0+^ or 
N"*"^. This means that low excitation species are relatively 
insignificant throughout the layer. Inspection of the data 
shows no evidence of broad wings underneath the low exci- 
tation lines at a level comparable to that observed in Ha, 
H/3 or [O III] . We have also explored whether faint wings 
would be present underneath the high excitation lines, as is 
predicted by the model. The limited S/N available for the 
weak [Ariv] lines, however, has made such a test inconclu- 
sive. 

The shapes of the profiles predicted by the TML models 
are either characterized by a broad profile with a peaked 
center ([O ill] , C iv, C ill] , [Ne ill] , Ha, H/3) or a narrow core 
without any significant wings ([Six], Mgii, [Nil], [On]). 
There are a few exceptions, however, such as the high exci- 
tation line of He 11 A4686, which presents a top hat shape. 
This arises as the line is only produced significantly at the 
onset of the thermal bump in Fig. [3] where the temperature 
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Fig. 8. Same notation as in Fig. [5] except for the 
TML model, which corresponds to a slower wind 
of T4)=3500kms^^ with A''2o=31. For comparison, a 
Lorentzian profile of 1800 km s~^ fwhm is overlayed. 

shoots up above that of the isothermal regior0 (see He 11 
emission fiux in Fig.|3]). For illustrative purposes we show 
the profile predicted for Hen in Fig. [51 It corresponds to 
an upper limit equivalent to 5 times the intensity predicted 
by the TML model with K,=3500kms^i and N°q=31. It 
is possible that mappings ic underpredicts the strength of 
this line, since it considers only pure recombination and 
does not include coUisional excitation. It is interesting to 
note that one cannot completely rule out the presence of 
a weak broad component (or flat pedestal) underneath the 
narrow Hell profile. We do not expect such a component 
to be due to a Wolf-Rayet feature, since supercluster A in 
NGC 2363 is much too young (< 1 Myr; Drissen et al. 2000) 
and, furthermore, no evidence has been found of broad N ill] 
or C IV, which are otherwise expected if Wolf-Rayet stars 
were present. 

4.4. Line profiles from starburst galaxies 

A different area where TML emission could play a 
role are broad profiles, with fwhm in the range of 
150-400 km s~^, observed in starburst galaxies, such as 
NGC 1569 (Westmoquette etal. 2007a, b, 2008) or M82 
(Westmoquette etal. 2007c). The requirements for wind 
velocities in these objects are modest in appearance, which 
implies a much lower wind velocity. However, the fwhm ob- 
served may nevertheless correspond to a significantly higher 
Vu, if the wind, rather than being fully isotropic, tended to 
project along the plane of the sky of the observer. For in- 
stance, in the case of near edge-on galaxies, the wind can be 
funneled above and below the galactic plane, resulting in a 

* Since the density decreases inversely to the temperature, the 
HellA4686 line emissivity becomes negligible towards the very 
high temperature end. 
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Fig. 9. Observed emission line profiles of He ii A4686 (black 
line) as a function of Doppler velocity with respect to the 
centroid of both lines. Other line transitions near He ii 
have been blanketed by dotting the corresponding line. We 
overlay a 3D profile from a TML model (gray line) with 
l/u)=3500kms~"'^ and N2Q—31, but with its intensity multi- 
pHed by a factor 5. 



narrower profile because 9i > 0, in the notation of Sect.| 
In contrast to NGC 2363, supernovae can play an important 
role in launching starburst winds. Another significant dif- 
ference is that it is the lower excitation lines like [N ii] A6583 
that are observed to be broad. This can only be accounted 
for in our model if a sufficiently low ionization parameter 
is assumed, which could ensure that the intermediate exci- 
tation lines spread over a significant fraction of the TML's 
thickness. A softer ionizing continuum due to high metallic- 
ities may also contribute to the prominence of intermediate 
excitation lines. In the case of M82, Westmoquette etal. 
(2007c) report a relatively low [O iii]/H/3 (5007/4861) ratio 
of ~ 0.32. New TML calculations adapted to the context 
of starburst galaxies would be required to verify whether 
the inferred ionization parameter is sufhciently low for the 
[N ll] lines to become broad as well. 

5. Conclusions 

Four mechanisms have been explored by Roy etal. (1992) 
to account for the extreme gas velocities observed in 
NGC 2363, namely stellar winds, Thomson scattering by 
hot electrons, supernova remnants, and a superbubble 
blowout. These authors reported significant problems with 
each of these mechanisms and considered them to be un- 
satisfactory. They concluded that the broad nebular gas is 
probably due to "very high velocity gas whose origin is, at 
present, unknown". In this paper, we explore the possibility 
that the broad emission originates from turbulent mixing 
at the interface between a hot cluster wind and more qui- 
escent photoionized condensations. In our model, all gas 



phases are in pressure equilibrium and, because of the high 
temperature of the wind, the total ionized gas mass con- 
tained in the wind is much lower than postulated by Roy 
etal. (1992), by a factor ^ Tyj/W^. Since the upper limits 
on Tu, set by the non-detection in the X-rays allows tem- 
perature values as high as 10* °K (Sect. 13. 3|) . the objection 
to the wind hypothesis expressed by Roy et al. can be lifted. 

The basic input parameters affecting the TML models 
have been given appropriate values as follows: we inferred 
a density of ~ 100 cm^^ for the nebular gas, using the [S ii] 
doublet. The metallicity that was assumed is 20% solar, 
as deduced by Luridiana etal. (1999). The SED that we 
adopted was calculated using the code LavalSB and corre- 
sponds to a newly formed star cluster 1 Myr old [the upper 
age limit derived by Drissen etal. (2000)]. Finally, the ion- 
ization parameter has been adjusted so that the calculated 
[Oiii]/H/? ratio became comparable to the observed value. 

In our model, the broad profile results from radial accel- 
eration in 3D of photoionized turbulent gas. About 10-20% 
of the acceleration that takes place occurs at temperatures 
approaching that of the wind, and this gas does not con- 
tribute to line emission. For this reason, the wind velocity 
assumed in the TML model must exceed the value inferred 
from the extent of the faint wings (i.e. ±3500 km s^^). To 
constrain the models, we adopted the spectral data set of 
Drissen et al. (2009) that was obtained at the Gemini ob- 
servatories. We report in Sect. 14. 21 on the results of mod- 
els with Vu, = 4300 km s^^, which can successfully fit the 
faint wings when the TML is ionization bounded (A''2o=33; 
Fig. [6]). The bulk of the profile is better fitted if we re- 
duced the TML thickness to N^^ ~ 10 (Fig.E]). This model, 
however, is matter-bounded and the [Oiii]/H/3 ratio from 
the TML would exceed that observed in the core of the 
lines as well as in the broad component. Our preferred sce- 
nario, therefore, is that the cluster wind undergoes radial 
acceleration, but the TML remains ionization-bounded at 
all radii. Our slower wind model of 3500 km s^^ and with 
iV2Q=31 is as successful in fitting the bulk of the broad 
profile (Fig. [8]). It is therefore our best model overall. The 
relative success of TMLs in fitting broad profiles and in 
predicting the absence of broad lines for [N ii] and [S ii] are 
significant arguments in favor of the models. The require- 
ment of wind velocities of order 4000 km s~^ in order to fit 
the very faint wings remain a stumbling block, however, un- 
less supercluster A, because of its extreme youth, is able to 
generate such physical conditions. Interestingly, our TML 
model results in a flat-top profile for He ii A4686, although 
this characteristics could not be confirmed with the cur- 
rent data. As a follow-up study, new calculations adapted 
to the broad component observed in starburst galaxies (e.g. 
Westmoquette et al. 2007a, b, c) should be carried out along 
the lines developed in Sect. 14. 41 

Although, in this paper, the accelerating mechanism is 
assumed to arise from a stellar wind, alternative mech- 
anisms could be considered that would lead to a near- 
constant radial acceleration and therefore generate profiles 
equivalent to those shown in this work. One possibility for 
instance might be radial acceleration due to radiation pres- 
sure acting on dusty photoionized gas plumes. NGC 2363 is 
clearly a fascinating object that remains a challenge in our 
understanding of extremely young nebulae. 
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